Screening of a macroion by multivalent ions: Correlation induced inversion of charge. 
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Screening of a strongly charged macroion by multivalent counterions is considered. It is shown 
that counterions form a strongly correlated liquid at the surface of the macroion. Cohesive energy of 
this liquid leads to additional attraction of counterions to the surface which is absent in conventional 
solutions of Poisson-Boltzmann equation. Away from the surface this attraction can be taken into 
account by a new boundary condition for the concentration of counterions near the surface. Poisson- 
Boltzmann equation is solved with this boundary condition for a charged flat surface, a cylinder 
and a sphere. In all three cases, screening is much stronger than in the conventional approach. 
At some critical exponentially small concentration of multivalent counterions in the solution they 
totally neutralize the surface charge at small distances from the surface. At larger concentrations 
they invert the sign the net macroion charge. Absolute value of the inverted charge density can be 
as large as 20% of that of the bare one. In particular, for a cylindrical macroion, it is shown that for 
screening by multivalent counterions predictions of the Onsager-Manning theory are quantitatively 
incorrect. The net charge density of the cylinder is smaller than their theory predicts and inverts 
sign with growing concentration of counterions. Moreover the condensation looses its universality 
and the net charge linear density depends on the bare one. 



PACS numbers: 77.84. Jd, 61.20.Qg, 61.25Hq 



I. INTRODUCTION 

Many objects with a much larger size than atomic 
are strongly charged in a water solution and are called 
macroions. One can think about a rigid polyelectrolyte 
which, in a water solution, dissociates into cylindrical 
macroion and monovalent small ions. DNA and actin 
are the best known examples of such biological polyelec- 
trolytes. Other important types of macroions are charged 
lipid membranes and charged spherical colloidal particles. 
Macroions are screened by smaller ions of the solution 
of both signs. A correct description of the screening of 
macroions is tremendously important for a calculation of 
properties of individual macroions, for example, the ef- 
fective charge or the bending rigidity. Screening also de- 
termines forces acting between macroions and both ther- 
modynamic and transport properties of their solutions. 

This paper examines the screening of a rigid macroion 
with a fixed and uniform distribution of charge on its 
surface. Three standard geometries are considered below 
— an infinite flat surface, an infinite cylinder and a small 
sphere — each uniformly charged with the surface density 
— a < 0. The standard approach for a description of such 
problems is the Poisson-Boltzmann equation (PBE) for 
the selfconsistent electrostatic potential tp(r) 
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ZjNpj exp 



(1) 



Here e is the charge of a proton, D ~ 80 is the dielec- 
tric constant of water, Zie is the charge of a small ion of 
sort i and N(n is their concentration at the point where 
■0 = 0. The number of papers using the analytical and 
numerical solutions of Eq. ([!]) is extremely largea. On 
the other hand there is an understanding that Eq. ([l]) 



neglects ion-ion correlations and is not exact. Devia- 
tions from the distribution of charge predicted by PBE 
were demonstrated numerically&u for the following prob- 
lem. Consider screening of a charged surface, x = 0, 
of a membrane or a film by a water solution occupying 
halfspace x > 0. Assume that there is only one sort of 
counterions with the charge Ze > and their concentra- 
tion N(x) = No exp(— Z eib/kBT) — > at x — > oo. In this 
case the solution of Eq. (yj) is very simple and has the 
Gouy-Chapman form 



N{x) = 
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(2) 



where A 
I = Z 2 Ib and Ib 



Ze/(2nla) is the Gouy-Chapman length, 
e 2 /(DksT) ~ 0.7 nm is the Bjerrum 
length. At large Z and a, the length A can become of the 
order of the size of the water molecule or even smaller. 
For example, at Z = 3 and a = 1.0 e/nm 2 A = 0.08 nm. 
This means that almost all ions are located in the first 
molecular layer at the surface or, in other words, they 
condense at the very surface of macroion. This raises 
questions about the role of their lateral correlations and 
the validity of the solution Eq. (|2|) . p, 

It was found by numerical methods^ that for a typical 
charge density a deviations from Eq. (|^) are not large 
for monovalent counterions, but they strongly increase 
with the charge of counterions Z. It was suggested in 
Ref. [|, ||, ^, that at Z > 2 repulsion between multi- 
valent counterions condensed at the surface is so strong 
that they form a two-dimensional strongly correlated liq- 
uid (SCL) in which the short order of counterions is sim- 
ilar to that of a Wigner crystal (WC) . This idea was used 
to demonstrate that two charged surfaces in the presence 
of multivalent counterions can attract each other at small 
distances. 
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A theory of the influence of SCL of multivalent coun- 
terions on the density of screening atmosphere of a 
macroion -has been suggested recently by Perel and 
ShklovskiiEl(PS). The PS's main idea is to treat sepa- 
rately two subsystems: two-dimensional SCL of multi- 
valent counterions at the very surface and their gas-like 
dilute phase at some distance to the right of the sur- 
face. In the SCL, PS explicitly take into account strong 
correlations using the energy of WC as a simple approx- 
imation for the free energy of SCL. On the other hand 
the gas-like phase is treated in the PBE approximation, 
while the effect of SCL is taken into account with the 
help of a new bounadary condition for PBE. 

In this article, PS theory is developed in several direc- 
tions. First, the phenomenon of charge inversion is stud- 
ied in greater detail and the inverted charge as function 
of the counterion concentration is found. In particular, 
the maximum possible value of the inverted charge is esti- 
mated. Second, PS approach is generalized to a spherical 
macroion. Third, I add the comprehensive discussion of 
approximations made in this theory. 

The next section starts with a review of the thermody- 
namic properties of a two-dimensional SCL, which were 
obtained by Monte-Carlo and other numerical methods. 
It is shown that for Z > 2 and typical charge density 
cr, SCL is characterized by a large negative chemical po- 
tential of ions. In other words, due to their lateral two- 
dimensional correlations ions are more strongly bound 
to the surface than in PBE approximation. This phe- 
nomenon can be understood as the attraction of a Z- 
valent counterion to its correlation hole in SCL. While 
PBE fails to describe these correlation effects at the sur- 
face, it works well at a distance from the surface, where 
the energy of attraction to the correlation hole is smaller 
than ksT and also N(x) is small enough that three- 
dimensional correlations are very weak. It is found below, 
that PBE becomes valid at x 3> 1/4 ^> A and that char- 
acteristic scales of PBE solution A 3> 1/4. Thus SCL 
together with the intermediate boundary layer, A < x < 
1/4, from the point of view of PBE, provide only a new 
boundary condition N(x) — N(0) at x = for the con- 
centration N(x) of multivalent counterions. It is derived 
in Sec. Ill from the condition of equilibrium of the gas- 
like phase with SCL. Due to the large negative chemical 
potential of SCL, the new boundary condition requires 
that N(Q) is exponentially small in the dilute phase. The 
section III also duscusses the structure of the intermedi- 
ate layer A < x < 1/4 between SCL and the dilute phase, 
where exponential decay of N(x) actually takes place. 

In Sees. IV- VI, PBE is solved with the boundary con- 
dition for N(0) for the standard problems of screening of 
a charged flat surface, a cylinder and a sphere for different 
salt compositions of the bulk solution. In planar geome- 
try and for the bulk concentration of Z : 1 salt N(oo) = 
I found that at x > 1/4 N(x) obeys Eq. (|J) in which A is 
replaced by exponentially large length A. At finite N(oo) 
consequences of correlation induced binding of counteri- 
ons to the surface become even stronger. Even at expo- 



nentially small N(oo) = N(0) the positive charge density 
of SCL completely compensates — cr. At N(oo) > N(0) 
it becomes even larger than cr, so that total charge of 
the surface becomes positive. This phenomenon is called 
charge-inversion. First, it was noticed in numerical calcu- 
lationsB. This paper presents a theory of charge inversion 
for the case of screening by small size ions, which is to- 
tally based on correlation effects. Recently a number of 
publications discussed similar pheBjfflnanon. fp r screening 
of macroions by charged polymersBlijOo'LJ. 

In the case of a cylinder, the conventional picture of 
nonlinear screening called the Onsager-Manning conden- 
sation should be strongly modified when dealing with 
multivalent ions. Consider a cylinder with a negative 
linear charge density — r\ and assume-that rj > ry c , where 
rjc = Ze/l. Onsager-Manxung theoryca, confirmed by the 
solution of PB equationtj, shows that such a strongly 
charged cylinder is partially screened by counterions re- 
siding at its surface, so that net (total) linear charge den- 
sity of the cylinder, ij* is equal to the negative universal 
value — rj c . The rest of the charge is screened at much 
larger distances according to linear Debye-Hiickel theory. 

The Onsager-Manning condensation does not take into 
account lateral correlations of counterions. In Sec.V the 
role of these correlations is considered and an analytical 
expression for 77* as a function of — r; and concentration 
A (00) of Z : 1 salt is derived (see Eq. (p4|)). It is shown 
that due to additional binding of multivalent counteri- 
ons provided by their SCL on the surface of cylinder, the 
absolute value of the negative net charge density, 77*, is 
smaller than in the Onsager-Manning theory. Moreover 
it strongly depends on a bare linear density, —77, so that 
attractive universality of the Onsager-Manning theory is 
destroyed. When concentration of counterions in the bulk 
N(oo) grows, the net charge density, rf changes sign from 
negative to positive at the point where N(oo) = N(0). 
Thus, charge inversion takes place for a cylindrical geom- 
etry, too. Positive 77* continues to grow with A(oo) until 
it reaches critical value, e/ls, for the Onsager-Manning 
condensation of monovalent negative ions. 

Finally, this paper studies screening of an uniformly 
charged small sphere with a negative charge —Q and ra- 
dius, a. For a strojiglvj-charged sphere the solution of 
PBE is well known EjOEj and is approximately valid for 
monovalent counterions. It shows that, in contrast with 
a charged plane or cylinder, a sphere has no condensed 
counterions, if AT (00) = 0. This happens because the po- 
tential energy of a counterion on the surface of totally 
ionized sphere, —QZe/Da, is finite. At N(oo) > N c , 
where N c oc exp(— QZls/ae) is an extremely small con- 
centration, a fraction of positive screening charge con- 
denses at the surface of the sphere and partially com- 
pensates its charge, so that the net charge of the sphere, 
Q* , changes in the range > Q* > —Q. In this regime, 
Q* does not depend on Q. This universality is similar to 
the one of the Onsager-Manning theory. The rest of the 
screening atmosphere can be described in Debye-Hiickel 
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approximation. When N(po) grows, Q* becomes smaller 
in absolute value but remains negative. 

It is shown in Sec. VI that in the case of screen- 
ing by multivalent counterions due to additional bind- 
ing by SCL, the net charge, Q* , behaves differently. 
When N(oo) < N c , all counterions are still lost. But 
at N(0) > AT(oo) > iV c a larger amount of counteri- 
ons condense at the sphere than PBE predicts. As a 
result at N(oo) — N(0) the net charge, Q* , vanishes. 
At N(oo) ^> N(0) positive Q* continues to grow and 
saturates at the value Q* ~ \JQZe. It should be empha- 
sized that when correlations are taken into account the 
above mentioned universality disappears and Q* becomes 
a function of the bare charge, Q. 

Note that the net linear charge density rf of a cylinder 
and the net charge, Q*, of a sphere are measurable quan- 
tities. In both cases they include only counterions which 
binding energy exceeds fc^T. These counterions move 
together with a cylinder or sphere, for example, in the 
electric field. Therefore, 77* and Q* can be studied in an 
electrophoresis experiment. In the case of charge inver- 
sion, a cylinder, sphere and any other macroion should 
drift in an anomalous direction. 

In Sec. VII approximations of this theory are dis- 
cussed. In the Conclusion several possible extensiones 
of this theory are mentioned. 



II. WIGNER CRYSTAL AND STRONGLY 
CORRELATED LIQUID 

It is shown below that for a > 1 e nm~ 2 and Z > 2, 
almost all charge of the plane is compensated by SCL of 
counterions at its surface, which has a two-dimensional 
concentration n = a/Ze. In this section I discuss ther- 
modynamic properties of this two-dimensional system. 
The minimum of the Coulomb energy of counterion mu- 
tial repulsion and their attraction to the background is 
provided by triangular close packed WC of counterions. 
Let us write energy per unit surface area of WC at T = 
as E — ne(n) where e(n) is the energy per ion. One 
can estimate e(n) as the interaction energy of an ion 
with its Wigner-Seitz cell of the background charge (a 
hexagon of background with charge —Ze and counterion 
in the center, which can be approximately viewed as a 
disc with radius R = (to) _1 ' 2 ). This estimate gives 
eCa,) ~ —Z 2 e 2 /DR. A more accurate expression for e(n) 



2 C 2 



Z 2 e 



e(n) = ~an 1/2 Z 2 e 2 D- 1 , 



(3) 



where a = 1.96. Eq. ([|) can be rewritten in units of the 
room temperature thermal energy, fc^T, as 



e(n) £ii -1.4 Z 3 / 2 (anm 2 /e) 1 /2fc B r. 



(4) 



The inverse dimensionless temperature of SCL is usually 
written in units 



RDk B T 



= 0.9 



Hn)\ 
k B T 



(5) 



For example, at a — 1.0 e/nm 2 and room temperature, 
Eq. (g) gives T = 1.2, 3.5, 6.4, 9.9 at Z = 1, 2, 3, 4. 
Thus, for multivalent counterions one deals with a low 
temperature situation. T is the large parameter of this 
paper theory. In its terms R/X ~ 2r 3> 1 and l/R ~ T 3> 
1. For example, at Z = 3 and a = 1.0 e/nm 2 lengths 
A, R and I are equal to 0.08, 1.0, 6.3 nm respectively. 
The small value of A means that almost all counterions 
are located in the first molecular layer at the surface and 
literally form a two-dimensional system. 

It is known, however, that, due to the smalLshear 
modulus, WC melts at a very low temperature^^ near 
r ~ 130. Nevertheless, the disappearance of the long 
range order only slightly changes thermodynamic prop- 
erties of the system. They are determined by the short 
range order, which, in the range 5 < T < 15 r iShQiild not 
be significantly different from that of the WCu'B'QtjLT-liiS| 
statement is confirmed by numerical calculational3l23'Eil 
of thermodynamic properties of the two-dimensional SCL 
of Coulomb particles on the neutralizing background or 
so-called one-component plasma. In the large range, 
0.05 < T < 5000, the excess internal energy of SCL 
per counterion (difference between internal energy and 
energy of the ideal gas with the same concentration) , 
e(n,T) = fcsT/(r), was fitted by the expressionllSl 



f(T) = -l.ir + 0.58r 1/4 - 0.26 



(6) 



with an error less than 8% (less than 2% in the range 
0.5 < r < 5000). The first term on the right side of 
Eq. (||) dominates at large T and leads to Eq. (||). The 
other two terms provide a relatively small correction to 
the energy of WC. It is equal to 11% at T — 5 and to 5% 
at T = 15. The reason for a such small correction is that 
short range order in SCL is similar to that of WC. For 
the free energy of unit area, F, one can write 



F = F(T = 0.05) + nk B T [ f(T')dr'/T', 

J0.05 



(7) 



so that for the chemical potential which is used below to 
describe the equilibrium of SCL with the gas-like phase 
one obtains 

fi(n, T) = -k B T\n(n w /n) + fj, s + /j, c (n, T), (8) 
fi c (n,T) = -fc s T(1.65r-2.6ir 1/4 + 0.26 lnT + 1.95). (9) 

Here fi c is contribution of correlations to the chemical 
potential. The high temperature fj,(T — 0.05) with suf- 
ficient accuracy is replaced by the chemical potential 
—kBT\n(n w /n) + ji s of an ideal two-dimensional solu- 
tion of ions in the surface layer of water with a two- 
dimensional concentration n w . The term /i s is the hy- 
dration free energy per ion at the surface and atpH. <C n w 
does not depend on the concentration of ions n c3. 
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The first term of Eq. (||) corresponds to the WC pic- 
ture. Indeed, one can find directly from Eq. (0) and 
Eq. (§) that 

d\ne(n)] 3 , . „, „ . 

mc = 1 , V n = r£ (n) = -1.65 rfc B T. 10 
an 2 

At large T, the chemical potential /i^C dominates in 
Eq. (||). The last three terms of /j c give 20% correction 
to the WC term at T = 5 and only 10% correction at 
r = 15. Thus, if necessary, at 5 < T < 15 one can 
use nwc as a first approximation. Below fi c is always 
calculated using the full Eq. (g). 

All necessary information about two-dimensional SCL 
has been presented. It is time now to study its equilib- 
rium with the rest of the screening atmosphere. 

III. A NEW BOUNDARY CONDITION FOR 
POISSON-BOLTZMANN EQUATION 

When an ion moves away from SCL it leaves behind 
its negatively charged correlation hole. If U(x) is the 
correlation energy of attraction to the hole, condition of 
equilibrium between SCL at x — and the gas-like phase 
at a distance x can be written as 

fi(n) + Zeip{0) = n(N) + Z&p(x) + U{x), (11) 

Here /x(n) is given by Eq. (|^), Zeip(x) is the counterion 
energy in the selfconsistent potential and 

fi(N) = -k B TlnN w /N + fib (12) 

is the chemical potential of the bulk gas-like phase, N w 
is the bulk concentration p£ water and fib is the is the 
bulk hydration free cncrgyEa, which does not depend on 
N . According to the terminology of Ref. |23| Eq. jll| ) 
means that the electrochemical potential of counterions 
is constant. 

It will be shown below that U(x) becomes less than 
k B T at x > 1/4. On the other hand, in many impor- 
tant cases the surface is screened so strongly that the 
selfconsistent potential changes by k B T only at exponen- 
tially large length, A, which is defined below. Therefore, 
the condition of equilibrium between SCL and the layer 
1/4 < x < A is 

M (n)=/z(A0> (13) 

Using Eq. (§) and Eq. @ and solving Eq. for N(x) 
one obtains that at 1/4 <§; x <C A concentration N(x) 
does not depend on x and equals 

^(OH-exp(-M^), (14) 
w \ k B l J 

where w — (n w /N w ) exp[(/ib — fj, s )/k B T)]. Below it is 
assumed for simplicity that /ij = /x s , i. e. surface and 
bulk hydration free energies are equal. In this case w is 



the length of the order of size of the water molecule (for 
estimates w — 0.3 nm is used below). 

The notation N(Q) reflects the fact that this value plays 
the role of a new boundary condition at x -C A for im- 
portant solutions of PBE which have large characteristic 
length A ^> 1/4. In such class of solutions N(Q) provides 
an universal description of the role of SCL. This paper 
deals only with this class of problems. Situations when 
one has to go beyond the universal boundary condition 
( |l4|) and start directly from Eq. ( |Tl|) will be studied in 
the next paper. 

Due to the large value of \/j, c (n,T)\, the concentration 
N(0) can be very small. For example, at a = 1.0 e/nm" 2 , 
and Z — 2,3 and 4 at which F =3.5, 6.4 and 9.9, and one 
gets, according to Eq. (g), \fi c (n,T)\/k B T =4.5, 8.8 and 
14.3 respectively. This gives for N(0) = 30 mM, 0.3 mM 
and 0.8 fiM for Z = 2, 3 and 4 (1 M = 6 10 26 m~ 3 ). It 
is clear now that |/i c (n, T)| plays the role similar to the 
work function for thermal emission, to the free energy of 
chemosorption or to the evaporation energy for the cases 
of equilibrium gas-liquid or gas-solid interfaces. The con- 
centration N(0) is similar to the density of the saturated 
vapor. 

Thus, correlation effects in SCL provide additional 
strong binding of counterions to the macroion surface. 
Qualitative arguments for such binding can be found in 
Ref. ||. We would like to stress that such binding does 
not happen at Z = 1. Indeed, at a = 1.0 e/nm -2 
one obtains from Eq. (p|) and Eq. (0) that r = 1.2 and 
fJ, c ,(n, T)/k B T = 1.3. Therefore, the boundary condition 
Eq. (|1J) does not produce nontrivial effects and standard 
solutions of PBE remain approximately valid. 

Below, I justify the role of the distance 1/4 and give 
an idea how N(x) evolves from n/A at x ~ A to N(0) 
at x — 1/4. Let us move one ion of the crystal along 
the x axis. As it is mentioned above, it leaves behind its 
correlation hole. In the range of distances A <C x <C R, 
the correlation hole is approximately a disc of the surface 
charge with radius R (the Wigner-Seitz cell) and the ion 
is attracted to the surface by its uniform electric field 
E = 2na/D. Therefore, if A were larger than w one 
would get N(x) = (n/A) exp(— x/X) at x <C R. In the 
cases of our interest A < w and N = n/w at x < w, 
while at w <C x <C R 

Ti 

N(x) = - exp(-x/X). (15) 
to 

At x 3> R the correlation hole radius grows and becomes 
of the order of x. Indeed SCL on the uniform background 
can be considered as a good conductor in the plane (y, z). 
It is known that a charge at a distance, x, from a metallic 
plane attracts an opposite charge into a disc with the ra- 
dius ~ x or, in other words, creates its point-like image 
on the other side of the plane at the distance 2x from 
the original charge. The same thing happens to SCL. 
The removed ion repels other ions of SCL and creates a 
correlation hole in the form of a negative disc with the 
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charge — Ze and the radius x. The correlation hole at- 
tracts the removed ion and decreases its potential energy 
by the Coulomb term 



U{x) = -ZV/4Da: 



(16) 



This effect provides the correction to the activation en- 
ergy of N(x) : 

* * tl ( \fi(n)\- Z 2 e 2 /ADx\ . . . . 

jV(x) = -exp^- ^ ^ J (x^R). (17) 

The similar "image" correction to the work function of 
a metal is well-known in the theory of thermal emission. 
The correction decreases with x and at x = I /A, becomes 
equal to fcgT, so that N(x) saturates at the value N(0). 
The dramatic difference between the exponential decay 
of Eqs. (15) and jl7| ) and the 1/x 2 law of Eq. (0) is obvi- 
ously related to the correlation effects neglected in PBE. 
Recall that it was assumed in the beginning of this paper 
that the charge of the surface is almost totally compen- 
sated by SCL. Exponential decay of N(x) with x con- 
firms this assumption and at T > 1 makes this theory 
self consistent. 

Consider now what happens with N(x) at distances 
i » i/4. At such distances, correlations of the removed 
ion with its correlation hole in SCL are not important 
and the correlation between ions of the gas phase are even 
weaker because N(x) is exponentially small. Therefore, 
one can return to PBE. In the next section solutions of 
PBE for the planar geometry for different concentration 
of salt are discussed. 



IV. PLANAR GEOMETRY 

The solution of PBE with the boundary condition ( |l"4| ) 
and N(oo) — is similar to Eq. (0): 



N(x) 



1 1 

2~7rl (A + x) 



0z>Z/4), 



(18) 



where the new renormalized Gouy-Chapman length, A, 
is exponentially large 



A = (2irlN(0)) 



-1/2 _ 



" exp (M^). (19) 



2iml 



2k B T 



For example, at a — 1.0 e/nm 2 , Eq. ( [l9| ) gives A ~ 
1.8, 12.3, 166 nm at Z = 2, 3, 4. These lengths should be 
compared with l/A — 0.7, 1.6, 2.8 nm respectively. We see 
that A > l/A for Z > 2. This justifies the use of Eq. (|l|) 
as the boundary condition for the large distance solution 
of PBE. 

Using Eq. (|l^) one finds that the total surface charge 
density located at distances x < l/A 



For a — 1.0 e/nm 2 one finds that cr* = 7 10 2 a at 
Z = 2, cr* = 7 10- 3 a at Z = 3 and a* = 4 10~ 4 rr at 
Z = 4. Corrections to fj, c (n,T) and N(0) related to such 
small a* can be, of course, neglected. 

One can compare these results with predictions of 
Eq. (||). Integrating Eq. (Bj from l/A to go one finds 
cr* = 2Ze/nl 2 , i.e. a* = 5 10" 2 e/nm" 2 at Z = 3 and 
a* = 2 10 -2 e/nm -2 at Z = A. These values of a* are 
much larger than Eq. (|20|). Thus, binding to the surface 
is strongly enhanced by correlation effects. 

Until now this paper has addressed the case of ex- 
tremely dilute solution, when N(oo) = 0. Consider the 
case of a finite concentration, N(oo), of a Z : 1 salt 
in the bulk of solution, or in other words, of a concen- 
tration N(oo) of Z-valent countcrions and concentration 
N-(oo) — ZN{po) of neutralizing ions with the charge 
— e. This adds the Debye-Hiickel screening radius 



r s = (4ttWV(oo)(1 + l/Z)) 



-1/2 



(21) 



to the problem. If N(oo) <C ^(0) the screening radius 
r s 3> A, and the fact that N(oo) is finite changes only 
the very tail of Eq. (18) making the decay of N(x) at 
x » r s exponential. At x <C r s , still N(x) ^S> ZN_{x) 
and all previous results are valid. However, when iV(oo) 
approaches N(0), the solution changes dramatically and 
a* vanishes. Indeed, when N(oo) — N(0) concentra- 
tion N(x) = N(oo)exp(-Zei}/k B T) = N(0) stays con- 
stant and potential ip(x) — at x > l/A. This means 
that the surface is completely neutralized at distances 
< x < l/A. 

If N(oo) N(0) negative charges dominate at x r s . 
Indeed in the PBE approach, 



N(x) = N(oo)exp(-Zetp/k B T), 
AT_(x) = N-{oo)exp(eip/k B T), 



(22) 
(23) 



and when concentration N(x) decreases with decreasing 
x, the electrostatic potential,-0(:r), grows and N-(x) in- 
creases. One can derive a boundary condition for N^(x) 
at x = from Eqs. (|2|) and (p3|). For this purpose, 
one should first express ip(0) through N(0) with help of 
Eq. (H), and then find 7V_(0) from Eq. (§3). This gives 



AT_(0) = ZN(oo)[N(oo)/N(0)} 1/z , 



(24) 



where N(0) is given by Eq. (|14[). Then the solution of 
PBE for iV_ (x) at x <C r s has the form similar to Eq. Jig) 



N_(x) 



1 



2ttI b (A_ +x) 2 



where 



A_ = (27r/ B ^_(0)) 



-1/2 



(25) 



(26) 



a* = -y/N(0)/(2nl B ) = —a(X/A). 



(20) 



and l B = e 2 / (Dk B T) is the Bjerrum length. To compen- 
sate the bulk negative charge the positive surface charge 
density of SCL becomes larger than a, so that the net 
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surface charge density a* becomes positive. Similarly to 
Eq. ©, it is 



'M(O) 



2-kI b 2ttI b A„ 



(27) 



This phenomenon is called charge inversion and is, of 
course, impossible in the framework of the standard PBE. 
Technically, charge inversion follows from the small value 
of N(0) in Eq. ([14]). Its physics is related to the strong 
binding of counterions at the charged surface due to the 
formation of SCL. Remarkably, when r » 1, this phe- 
nomenon happens under the influence of a very small 
concentration of salt. 

According to Eqs. (|^) and ( gjj) the net density a* 
continues to grow with N(oo) at iV(oo) 3> N(0). It is in- 
teresting to study how far it can grow and how strong the 
charge inversion can be. The use of PBE with the bound- 
ary conditions Eqs. (|l4|) and (|2|) is valid if A_ > I /4. 

To estimate the maximum value of <r*, which can be 
reached within the range of validity of this theory one 
can substitute A_ = 1/4 into Eq. (|27]). This gives 



2e 



2R 2 

r zW B 



(28) 



For example, at Z = 3 and a = 1.0 e/nm -2 one ob- 
tains a* ~ 0.15cr. To find a* as a function of N(oo) 
in the whole range where A_ ^S> 1/4, one should solve 
Eq. ( p7| ) selfconsistently substituting n = (a + a*)/Ze 
into Eqs. (0), (|J) and @). 

One can show that at R < A_ < I /4, when our theory 
based on the universal boundary condition Eq. ( |l4| ) is 
not valid, a* continues to grow. If A_ becomes smaller 
than the radius of a Wigner-Seitz cell R, negative ions 
screen each counterion separately. The effective charge of 
counterions becomes smaller than Z. This weakens their 
lateral interactions and makes N(0) larger. Therefore, a* 
starts to decrease. The maximum value of a* is reached 
at A_ ~ R and is close to e/ (2itIbR). For Z = 3 and 
a = 1.0 e/nm~ 2 this gives a* ~ 0.24a. 

I will not try here to make the above estimates of the 
maximum a* more accurate because of sensitivity of this 
estimate to the ion size due to the following reason. It 
was assumed above that when cr^^ is reached, all salt 
molecules are still fully dissociated in water, so that the 
concentration, N(oo) of ions with charge Z is equal to 
the concentration of the salt, N s (oo). In reality, at very 
large N s (oo), the concentration of fully ionized counteri- 
ons, N(x), saturates at the level 



N max (oo) ~ b- 3 eM~Ze 2 /bDk B T), 



(29) 



where —Ze 2 /bD is the Coulomb interaction energy of 
positive Z-valent ion with negative monovalent one at 
the minimum distance between them, b. In this case, the 
majority of counterions keep a negative ion. One can 
refer to such a complex as (Z — l)-ion. The transition 



to such regime happens when the concentration of salt, 
N s (oo), reaches N max (oo). 

Substituting N rnax (oo) into Eq. j25| ) and then Eq. j25|) 
into Eq. (|27j ) one finds that at Z = 3, a = 1.0 e/nm 
and b > 0.4 nm this limitation of dissociation is not im- 
portant. For smaller b charge density, a*, saturates at 
the value N max (oo) which is smaller than Eq. (p9|), and 
stays at this level until the concentration of (Z — l)-ions 
becomes so large that they replace fully ionized ions at 
the surface. This leads to the drop of a* . 

Note once more that dramatic changes of the screen- 
ing atmosphere described above do not happen at Z = 1 
when r ~ 1 and |/i c (n, T)\/ksT ~ 1. The standard 
Gouy-Chapman solution of the PBE, Eq. (||), remains 
valid in this case. 



V. SCREENING OF UNIFORMLY CHARGED 
CYLINDER 

Consider screening of an infinite rigid cylinder with 
a radius a, a negative surface charge density — a or, 
in other words, with a negative linear charge density 
— r\ = — 2ixaa. Assume that a is large enough so that the 
surface of the cylinder is covered by a two-dimensional 
SCL with R < 2na and with T ^> 1. Such a cylinder can 
be a first order approximation for the double helix DNA, 
where a = 1 nm, rj = 5.9 e/nm, a — 0.94 e/nm 2 , and 
for Z = 3 radius of the Wigner-Seitz cell R ~ 1 nm and 
I = 6.3 nm. 

A screening atmosphere of a cylinder is described by 
the concentration N(r), where r is the distanpepfrom the 
cylinder axis. The solution of PBE is knownEJld to cop. 
firm the main features of the famous Onsager-Manningtj 
picture of the counterion condensation. This solution 
depends on the relation between r\ and r\ c = Ze/l = 
ksT/eZD. For a weakly charged cylinder with rj <C t] c , 
the screening is linear and can be described by the Debye- 
Hiickel approximation. For rj > r/ c , screening becomes 
nonlinear and most of the screening charge, — (?) — rj c ), is 
located at the cylinder surface, while at N(oo) = the 
rest of the screening charge, rj c is spread in the bulk of 
the solution. This means that at large distances, the net 
charge density of the cylinder, rj*, equals —ij c and does 
not depend on rj. Note that this is different from the pla- 
nar geometry where all the charge is bound to the surface, 
so that far enough from the surface, the net surface den- 
sity vanishes (the finite a* was defined at x < 1/4). At 
a finite iV(oo), the charge density rj* is screened only at 
linear screening radius r s . 

Here I deal with strongly charged cylinder, rj 3> i] c . 
It easy to check that this inequality follows from our as- 
sumptions R < 2iva and r> 1. It is also fulfilled for the 
case of DNA for which rj/rj c ~ 4Z. The goal in this case 
is to verify whether in the case of multivalent counteri. 
ons elegant statements of the Onsager-Manning theorytil 
should be changed due to SCL at the surface of a strongly 
charged cylinder. 



G 



As in the previous section, the boundary condition 
Eq. ([lj) is used below to allow for additional binding 
of counterions by SCL. One can introduce a radius rr, 
at which energy of interaction between a counterion and 
its correlation hole, U(r), becomes of the order of /j^T, 
so that the boundary condition N(r) = N(0) can be 
used. For a cylindrical geometry, ry, strictly speaking, 
differs from its analog for a planar problem 1/4. Indeed, 
at r ^> R energy U(r) can be calculated as energy of 
attraction the charge Ze to an infinite metallic wire with 
the radius a: 



U(r) = -Z 2 e 2 /4D(r-a) 
U(r) = -TrZ 2 e 2 /4Dr\n(r/a) 



(R < r - a < a) 
(r > a). 



(30) 



One can find tt from the equation |{7(rr)| 

(1/4 < a), 
(7/4 > a). 



k B T: 



r T = 
r T = 



■1/4 
kI 



41n(Z/4a) 



(31) 



At distances < r < r s the electrostatic potential of the 
linear charge density tf is not screened and the boundary 
condition of Eq. ( |l4| ) can be used to write 



N(r) 



N(0) exp 
A(0) exp 



Ze[ip(r) - ip(r T )} 



2rj* 



\n{r/rT 



(32) 



At r = r s concentration N(r, 
this equation for rf is 



N{oo) 

ln[A(0)/A(oo)] 
]n(r s /rT) 2 



The solution of 



(33) 



According to Eq. (pi]), at a not very largej/4a one can use 
estimate, tt ~ 1/4. Substituting Eq. pi] ) into Eq. (33), 
one arrives at 



V 



-Vc 



ln[A(0)/A(oo)] 
ln[4/(7rAT(oo)/ 3 )]' 



(34) 



It is clear from Eq. ( p4| ) that if two logarithms are close 
to each other, i. e. if 



In n / r > 1 

AT(oc) 



(35) 



the Onsager-Manning theory is approximately correct 
and rj* approaches — r\ c . If T ~ 1 and fj, c (n,T) ~ k B T 
concentration N(0) ~ n/w is large and inequality ( |3^ ) is 
fulfilled at any reasonable N(oo). Thus for typical charge 
density a and Z = 1 the Onsager-Manning result is red- 
erived. However, for screening by multivalent ions r>l, 
in Eq. ([u]) fi c (n,T) > k B T and concentration N(0) is 
exponentially small so that values of N(oo) at which rf 
is close to —rj c are extremely small. For example, to get 
rf = -0.75r] c one needs N(oo) = N . 75 = 0.02N(0) 2 l 3 . 



At a = 1.0 e/nm~ 2 and Z 
that 2V(0) = 1.7 1 23 m 



-3 _ 



= 3 it is shown above 
0.3 mM and therefore 



N , 75 = 2 10 2U m- d = 0.3 iiM. Switching to Z = 4, one 
has N(0) = 5 1 20 m" 3 = 0.8 which results in an 
unrealistically small Afo.75 = 2.5 10 15 m~ 3 . 

On the other hand, in disagreement with the Onsager- 
Manning theory one obtains from Eq. (Q) that <^ rj c 
when a concentration N(oo) of the salt is still exponen- 
tially small, namely N(0) 2 l 3 < N(oo) < N(Q). More- 
over, according to Eq. ( p4| ) 77* vanishes at N(oo) = N(0). 
This result is easy to understand without calculations. 
Indeed, in this case N(r) = N(oo) exp(—Zeip(r)/kBT) = 
N(0) stays constant and ip(r) = at all r > 1/4, so that 
all of the charge of the polyelectrolyte is compensated 
inside cylinder with r = 1/4. 

The difference from the Onsager-Manning theory be- 
comes even more apparent at N{po) > N(0) when the 
density 77* becomes positive. Note that this charge in- 
version takes place still at exponentially small AT (00). A 
positive rf continues to grow with N(po) until it reaches 
critical density 



e/h 



(36) 



and the standard Onsager-Manning condensation of 
monovalent negative ions starts. According to Eq. ( |3^ ) 
this happens at A (00) = N satl where 



N sa t ~ l~ 3 (N(0)l 3 ) 1 ^ z+1 \ 



(37) 



At N(oo) > N sat charge density 77* remains fixed at the 
level e/ls- Condensed negative ions eventually screen 
lateral interaction of counterions in SCL, \^ c (n,T)\ de- 
creases and 77* drops. Comparing Eq. ( |37j ) with Eq. ( |29| ) 
for the maximum concentration, A^ max (oo), of fully dis- 
sociated Z : 1 salt one sees that they are quite close, if 
ion size b is not too small. For a very small distance of 
the closest approach b, the growth of 77* is limited earlier 
by the condition N(oo) < N max (oo). 

To summarize, the net charge 77* as a function of salt 
concentration N(oo) is given by Eq. (|34|). It changes in 
the range ejls > rf > —r\ c when A(oo) grows. Strictly 
speaking, to quantitatively describe 77 as a function of 
iV(oo), one should use the self consistent concentration, 
n = (77 + 77*)/27ra, in Eq. @) for N(0). 

Finally, it should be emphasized that this result 
does not demonstrate the attractive universality of the 
Onsager-Manning theory. According to Eq. (|34|), 77* de- 
pends on 77 through fi(n) in Eq. ([li]). Thus, for the 
screening by multivalent ions at T > 1, and at any rea- 
sonable N(oo) the predictions of Refs. |l4|, [[5] are quali- 
tatively incorrect. 

Return now to the case of a small concentration of 
a Z : 1 salt, N(oo) <C A(0), and consider what hap- 
pens to 77* when a 1:1 salt with a larger concentra- 
tion, Ai ^> N(oo), is added to the solution. This is 
a realistic experimental situation. Such a problem can 
be solved with the help of Eq. ( ^3|), if one substitutes 
r s = (8-kIbNi)- 1 / 2 instead of Eq. (0). The result is 
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V = ~Vc 



ln[jV(0)/jV(oo)] 
ln[2^ 2 /(7riViZ 3 )]' 



(38) 



Eq. (||) shows how at a given N(oo) < N(0), the ab- 
solute value of the net negative charge density increases 
with N\. At large enough Ni, monovalent counterions 
replace counterions with charge Ze at the surface of the 
macroion. This replacement happens when correspond- 
ing change of free energy vanishes, i.e. at 

Ho{n,T) - k B TMn w /n) - Zk B T In^/JVi) = 
Z\i x - Zk B T\n{n w /Zn) - k B T\n(N w /N(oo)). (39) 

Here the left side is the free energy of a Z-valent coun- 
terion at the surface and Z monovalent ones in the bulk, 
while the right side is the free energy of Z monovalent 
ions at the surface and a Z-valent ion in the bulk. u\ is 
the correlation part of the chemical potential of a mono- 
valent ion. One can neglect Zfii in comparison with 
/i c (n,T) because, as it was mentioned above, /ii is nu- 
merically small and the latter quantity is proportional to 
Z 3 / 2 . Solving Eq. <M) for Ni one finds 



Ni 



n ( N(oo) \ 1/Z 
w V N(0) 



(40) 



Substituting this Ni into Eq. (p8j), one sees that at the 
moment of replacement rf ~ —rj c Z = —k B T/e, provid- 
ing a natural crossover to the case of screening by exclu- 
sively monovalent counterions. 



VI. SCREENING OF AN UNIFORMLY 
CHARGED SPHERE 

Consider application of this theory to a sphere with 
a small radius a = 2 — 5 nm and with a charge — Q 
screened by Z : 1 salt with concentration N(oo) in the 
bulk. At large enough surface charge density a = ZeriQ = 
~Q/4na 2 the surface is covered by SCL of Z-valent coun- 
terions. The goal is to find a concentration n of this SCL 
and the net charge of the sphere, 



Q* = 4ira 2 nZe-Q, 



(41) 



as a function of a, Q, Z, T and N(oo). In the case of a 
sphere, screening atmosphere is characterized by concen- 
tration of Z-valent ions, N(r), as a function of distance 
r from the sphere center. For simplicity, assume that 
a > 1/4 so that the boundary condition Eq. ( |l4| ) is valid 
at r = a + 1/4 where the curvature of the sphere can be 
neglected. At distances from the surface r — a <C r s one 
can neglect screening and find the concentration, N(r), 
similarly to Eq. (B2) as 



N(r) = N(0) exp 
~ N(0) exp 



Ze[ip{r) - ip{a + 1/4)] 
k B T 
Q*Ze(a- 1 - r- 1 ) 
Dk B T 



(42) 



At the distance r — a + r s where linear Debye-Hiickel 
theory starts to work, one has 

N(a + r s ) ~ AT(oo). (43) 

Solving this equation in the case r 8 3> a one obtains 



<i< , ^V(O) a ( n/w | \[i c (n,T)\ 



In 



l B Z N(oo) l B Z V ^V(oo) 



k B T 



(44) 



This equation is similar to Eq. ( p4[) . It is to be solved for 
Q*, n and |^ c (n,T)| together with Eq. @ and Eq. (|) 
(or its low temperature version Eq. (|10|)). 

In the case of monovalent counterions, when r ~ 1, 
fi c (n,T)/k B T ~ 1, so that correlations do not play 
any role in Eq. (B) and solution does not differ from 
the solution of PBE. In this case, the concentration 
iV(0) ~ n/w ~ uq/w is much larger then any reasonable 
N(oo}_ so that ln(AT(0)/AT(oo)) > and Q* < 0. Thus, 
Eq. ( J44| ) describes the partial compensation of charge — Q 
by positive charge Q + Q* = 4ita? Zen of counterions con- 
densed at the very surface of the sphere. The rest of the 
screening charge, —Q* is situated at the distance r s from 
the sphere in Debye-Hiickel atmosphere. It is clear now 
that nonlinear screening of a sphere is similar to Qaa ager .- 
Manning condensation in the case of a cylinder li3'EZlli3. 
In both cases there are two separate groups of counteri- 
ons: condensed and free. Moreover, for a sphere there 
is a similar universality of the net charge Q* . Indeed, 
when N(oo) Crio/w the dependence of Q* on Q is neg- 
ligible. (One can evaluate this dependence substituting 
no for n in Eq. (fSj).) The only qualitative difference be- 
tween the screening of a sphere and a cylinder is that at 
unrealistically small A(oo) < N c where 



N r = 



1 



4na 2 w 



exp 



Ql B Z 



ae 



(45) 



the last counterion leaves the surfaceE! and Q* = —Q. 

On the other hand, in the case of screening by mul- 
tivalent ions, correlations significantly change the above 
described mean field dependence of Q* on N(oo). These 
changes start, however, only at N(oo) > N c , because 
without condensed counterions correlations can not play 
any role. At N(oo) ^S> N c , when n grows and be- 
comes comparable with n$ one obtains that T 3> 1, 
fi c (n,T) > k B T and, according to Eq. (Q), N(0) is ex- 
ponentially small. Therefore, it follows from Eq. (|44|) 
that the negative net charge, Q*, grows (decreases in the 
absolute value) faster than it does in the case of mono- 
valent. Eventually, Q* vanishes at exponentially small 
N(oo) = N(0). 

At 7V(oo) > N(0) the net charge, Q*, becomes positive 
and continues to grow. As in other geometries this charge 
inversion happens because of strong binding of counteri- 
ons by SCL. At large enough N(oo) one can neglect the 
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first term in parenthesis of Eg. (f44|). Then using a low 
temperature expression Eq. ( |l0[) for n c (n) one find that 
Q* saturates at the positive value 

Q* max = p^/QZ~e, (46) 

where (3 — 3a/4v / 7r ~ 0.84. For example, at Z = 3 for 
a sphere with bare charge — 50e (Q = 50e) one arrives 
at the net charge Q* max — +10e. Eq. (^6|) is remarkably 
simple and universal: Q* max does not depend on a sphere 
radius a. 

Eq. (f46| ) for Q* nax is valid until one of the two follow- 
ing events happens. First, a concentration N(oo) can 
reach maximum concentration N max (oo) of fully dissoci- 
ated Z : 1 salt (see Eq. (l37|)). Secondly, a condensation of 
monovalent negative ions on the positive sphere can start. 
Condensed counterions eventually screen the lateral in- 
teraction of counterions in SCL and effectively change 
their charge from Z to Z — 1. As a result |/i c (n,T)| 
decreases and q starts to drop. Similarly to Eq. (pTu 
condensation starts at N(oo) oc exp(— /i c /(Z + ljfcsT). 
This concentrations is close to N max (oo), if ion size, b, 
is not too small. For a very small b condition iV(oo) <C 
N max (oo) is more restrictive. For a small sphere with 
radius a = 2 — 5 nm both restrictions start to work while 
r s 3> a. Therefore, I do not consider here the case r s < a. 
For larger spheres, r s can become smaller than the above 
mentioned other limits on N(oo) start to work. In this 
case, the sphere effectively works as a flat surface and 
one can use results of Sec. IV. 



VII. DISCUSSION OF APPROXIMATIONS 

In this section, approximations used in this paper are 
discussed. First, it was assumed that charges at the sur- 
face of a macroion are fixed and can not move. In the 
case of a solid or glassy surface, for example, colloidal 
particles and rigid polyelectrolytes, such as double helix 
DNA, this approximation seems to work well. On the 
other hand, for charged lipid membranes it can be vio- 
lated. If the surface charges are mobile they can accu- 
mulate near Z-valent counterions, forming short dipolcs 
directed perpendicular to the surface. These dipoles in- 
teract weakly with each other so that the energy of their 
lateral correlations is smaller than in SCL on the uniform 
background. On the other hand, such concentration of 
the surface charge under the counterion by itself creates 
an additional binding of counterions to the surface. As a 
result the negative chemical potential, fi{n,T), becomes 
larger in the absolute value and the boundary concentra- 
tion N(0) becomes smaller. The theory of this paper ex- 
presses everything through N(0). Therefore, the unusual 
effects of complete compensation of macroion charge and 
of the charge inversion become stronger. 

The second approximation made above is the assump- 
tion that fixed charge is uniformly distributed at the sur- 
face. Localized charges are actually discrete. Therefore 



it makes sense to discus whether -e charges, for example, 
randomly distributed on the surface, work as an uniform 
background. In the limit Z> 1, the repulsion between 
Z-valent counterions is much stronger than their pinning 
by the surface charges, so that the concept of an uniform 
background works exactly. At Z > 3 the uniform back- 
ground is still a good approximation for realistic values of 
the radius of closest approach, b s , of counterions and dis- 
crete surface negative charges. On the other hand, this 
approximation can fail at Z = 1 because all counterions 
and discrete negative surface charges have a tendency to 
form neutral Bjerrum pairs instead of SCL if b s -C R. In 
this case, N(0) can be small even for Z = 1. When neg- 
ative charges are clustered, for example, form compact 
triplets, even at Z = 3 interaction with such a cluster 
can be as important as interaction with the neighboring 
counterions. Each counterion temds to neutralize one 
cluster forming a neutral dipole. Again, this leads to 
stronger binding to the surface and smaller N(0). 

This discussion naturally leads us to the third ap- 
proximation used above in the calculation of the chem- 
ical potential of SCL. Eq. (^|) was obtained for point- 
like counterions. Actually counterions have a finite size 
and one wonders how this affects these results. Our re- 
sults, of course, make sense only if the counterion ra- 
dius, b c is smaller than the radius of a Wigner-Seitz cell, 
_R, so that counterions occupy a small fraction of the 
surface. In other words, the idea of SCL or WC works 
only when objects of large charge density arrange them- 
selves on a background with much smaller charge den- 
sity. For treevalent ions on the surface with the charge 
density a — 1.0 e/nm" 2 the radius of Wigner-Seitz cell 
R = 1 nm, so that for a counterion with b c = 0.5 nm this 
condition is easily satisfied. Positive corrections to the 
energy per ion of WC are proportional (b c /R) 2 and ap- 
pear due to the fact that the charge finite size counterion 
can not be situated exactly in the potential minimum 
created by its nearest neighbors. 

Finally, all estimates in this paper are based on the use 
of dielectric constant of water D ~ 80. For the lateral 
interactions of counterions near the surface of an organic 
material with a low dielectric constant, the effective di- 
electric constant D can be substantially smaller. (In a 
macroscopic approach it is close to D/2). As a result, 
absolute values of |/i(n,T)| can grow significantly and 
N(0) may become even smaller. 



VIII. CONCLUSION 

In conclusion, the role of strong lateral correlations 
of Z-valent counterions condensed on a charged surface 
is studied. It is argued that a strongly correlated liq- 
uid (SCL), or in other words, a two-dimensional one- 
component plasma is a good model for these correlations. 
It is shown that, due to the additional binding provided 
by SCL, the concentration of counterions close to SCL 
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is exponentially small (see Eq. (p_4|)). This concentration 
depends only on Z and the surface charge density of the 
macroion a and serves as a boundary condition for the 
Poisson-Boltzmann equation (PBE), which is still valid 
far from the surface. PBE is solved with the boundary 
condition (|lj) for all three standard geometries. For a 
charged cylinder, it is shown that in the presence of SCL 
the Onsager-Manning condensation is strongly modified. 
The increasing bulk concentration of Z-valent counteri- 
ons, iV(oo), makes the net negative charge of the cylinder 
smaller than in the Onsager-Manning theory drives it 
through zero and makes it positive. Similar changes are 
predicted for a charged sphere with charge —Q. In this 
case, charge inversion can result in positive a net charge 
0.84^/QZe. All these phenomena happen while the con- 
centration of Z : 1 salt is still exponentially small. Tech- 



nically they follow from the boundary condition (14), 
which in turn is a result of a strong correlations of coun- 
terions of the surface layer. 

This theory can be applied to variety of other prob- 
lems. First, one can study more complicated solutions 
where a substantial concentration of 1 : 1 salt is added 
to Z : 1 salt. We gave only one example of such a prob- 
lem in the end of Sec. V. Second, this theory should be 
extended to a finite concentration of macroions. In this 
case problems of a global instability of such a solution 
should be addressed, too. Third, one can use a similar 
theory, for counterions of a larger size and nonspherical 
shapeEl, provided they have larger charge density than 
the macroion's surface. For example, a positively charged 
rigid flat surface can be screened by a_splution of a rod 
like polymer such as double helix DNAE£l. If projected to 
the plane, the negative surface density of DNA is by ab- 
solute value larger than the charge density of the plane, 
DNA rods form strongly correlated nematic liquid, which 
provides strong binding of DNA to the surface. A net 
charge density of the plane can experience correlation in- 
duced charge inversion. These and other problems will 
be addressed in future publications. 
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